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ABSTRACT 

The traditional view that Lya emission and dust should be mutually exclusive has been questioned more and more often; most notably, 
the observations of Lya emission from ULIRGs seem to counter this view. In this paper we seek to address the reverse question. How 
large a fraction of Lya- selected galaxies are ULIRGs? Using two samples of 24/25 Lya emitting galaxies at z = 0.3/2.3, we perform 
this test, including results at z = 3.1, and find that, whereas the ULIRG fraction at z = 3.1 is very small, it systematically increases 
towards lower redshifts. There is a hint that this evolution may be quite sudden and that it happens around a redshift of z ~ 2.5. After 
measuring the infrared luminosities of the Lya emitters, we find that they are in the normal to ULIRG range in the lower redshift 
sample, while the higher redshift galaxies all have luminosities in the ULIRG category. The Lya escape fractions for these infrared 
bright galaxies are in the range 1 - 100% in the z = 0.3 galaxies, but are very low in the z = 2.3 galaxies, 0.4% on average. The 
unobscured star formation rates are very high, ranging from 500 to more than 5000 M Q yr~', and the dust attenuation derived are in 
the range 0.0 < A v < 3.5. 
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1. Introduction 

Galaxies at high redshift come in all kinds of flavours. 
Depending on the search criteria, the sources found may be 
dusty or dust-free, more or less massive, star forming or qui- 
escent. Some of the most impressive beasts of the high redshift 
zoo are the sub-mm and ultra-luminous infrared galaxies: SMGs 
(Ivison et al. 1998, 2000, Blain et al. 2002, Chapman et al. 
2005) and ULIRGs (Sanders & Mirabel 1996). These galaxies 
show star formation rates in the hundreds or thousands of solar 
masses per year, with most of their light re-processed into the 
infrared or sub-mm wavelengths due to large amounts of dust. 
On the other end of the scale, Lya emitters are found (LAEs; 
M0ller & Warren 1998, Gronwall et al. 2007, Nilsson et al. 2007, 
Finkelstein et al. 2009, Ouchi et al. 2008, Grove et al. 2009) 
which have hitherto been generally considered blue and dust- 
free galaxies with moderate star formation rates of a few to a 
few times ten solar masses per year. 

Only a few publications have reported finding red and/or 
dusty Lya emitters (Stiavelli et al. 2001, Colbert et al. 2006, 
Lai et al. 2007, Finkelstein et al. 2009, Nilsson et al. 2009). 
Considering the traditional view that even small amounts of dust 
would quench any Lya emission (for a discussion, see Pritchet 
1994), the detection of the Lya line in dusty galaxies would seem 
surprising, but even more so the detection of Lya in a sample of 
sub-mm galaxies (Chapman et al. 2003, 2005). Likely explana- 
tions for the existence of dusty galaxies with Lya emission are 
either special geometrical alignments or strong local variations 
in the dust-to-gas ratios in a clumpy medium (Neufeld 1991, 
Hansen & Oh 2006). Observational evidence that Lya emission 
is commonly seen among sub-mm galaxies does not, however, 
constitute evidence for the opposite, since sub-mm galaxies are 



rare and because there are few ways other than Lya to deter- 
mine the redshift of very red galaxies. It follows that the over- 
lap reported so far could well be an observational selection bias. 
The fundamental questions that remain are therefore whether i) 
sub-mm galaxies and ULIRGs are common among Lya emit- 
ters and ii) whether the ULIRG fraction evolves with redshift. 
This Letter aims to address those two questions. We here present 
the infrared properties of two samples of galaxies found through 
their Lya emission at z — 0.3 and 2.3. We show that they have in- 
frared fluxes well into the ULIRG regime. This sample is unique, 
bridging the gap between low-mass star forming and large star- 
bursting galaxies. 

Throughout this paper, we assume a cosmology with Hq = 
72 km s" 1 Mpc- 1 , Q m = 0.3 and Q A = 0.7. 

2. Data 

The Lya candidate samples studied here come from Deharveng 
et al. (2008, z = 0.3, hereafter z03) and Nilsson et al. (2009, 
z = 2.3, hereafter z23). The z03 sample consists of a total of 
3 1 Lya emitting galaxies at 0. 1 < z < 0.4 in the ECDF-S and 
ELAIS-S fields, found in a survey volume of ~ 5 x 10 5 Mpc 3 . Of 
the 31, one is observed to be an AGN by Cowie et al. (2009) 
and 11 are not included in their classification. At z — 2.3, 
the sample includes 187 Lya emitting candidates at redshifts 
2.21 < z < 2.31, and are spread over a ~ 0.2 deg 2 area in the 
central COSMOS field (the survey volume is ~ 3 x 10 5 Mpc 3 ). 
Of the 187 candidates, 27 are considered AGN, based on detec- 
tions in public Chandra and/or XMM X-ray images. The ECDF- 
S, ELAIS-S and COSMOS fields are covered by the SWIRE 
(Lonsdale et al. 2003) and S-COSMOS (Sanders et al. 2007) 
Spitzer surveys. The photometry of all the z23 candidates, in all 
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Fig. 1. Histogram of infrared luminosities of all IR detected 
sources as derived with Eq. Q] At both redshifts, the hatched 
area marks genuinely certain LAE galaxies, while the open ar- 
eas are AGN or unclassified objects. Solid vertical lines mark the 
LIRG/ULIRG definitions. 

Spitzer bands from 3.6 - 24pm, will be presented in a forthcom- 
ing publication (Nilsson et al., in prep.). Here we focus on the 
results found for the photometry in the 8pm and 24pm (MIPS) 
bands for the z03 and z23 samples, respectively, corresponding 
to ~ 7/im in the restframe of both samples. 

For the z03 emitters, aperture photometry centred on the Lya 
coordinates was performed. The aperture radius was 4.5", se- 
lected to be consistent with the z23 results. The limiting sensi- 
tivity is a few times ten pjy. Of the 3 1 candidates, 24 are detected 
at > 3<r, of which Cowie et al. (2009) classified one as an AGN 
and 17 as galaxies. 

In COSMOS, both a deep and a wide survey has been pub- 
lished in the MIPS band. The deep survey covers ~ 33% of the 
field surveyed for Lya, and the sensitivities reached in the two 
surveys were 71 and 150 pjy, respectively (5cr). To find coun- 
terparts to the Lya emitters, the public catalogues were searched 
within 4 pixels (4.8", i.e. 0.8 x FWHM of the MIPS PSF) radii 
of each source and 25 counterparts were found. Of these, 16 are 
also X-ray detected, and are considered to be AGN based on 
their R band-to-X-ray flux ratios. The fluxes of the galaxies in 
the 8/24 pm bands are found in Table Q] 

3. Infrared properties of Lya emitters 

At z = 2.25, the MIPS 24/im band corresponds to restframe 5.9 - 
8.8 pm. Correspondingly, the 8pm IRAC band covers 4-8.8 pm 
for the z03 emitters. To convert this mid-infrared luminosity to 
the total infrared luminosity, we use the conversion of Chary & 
Elbaz (2001): 

L IR [8 - 1000pm] = 4.37^ x 10" 6 L^. (1) 

The values derived can be found in Table Q] and a histogram of 
the total infrared luminosities is shown in Fig. Q] At each red- 
shift, the infrared luminosities for the full samples (25 and 24, 
respectively), as well as the subsamples with certain identifica- 
tions, are shown. It is seen that the flux limits at the two red- 
shifts cause the overlap between the high and z03 sample to be 
very small. All of the z23 sources are consistent with ULIRG 
luminosities (L/« > 1O 12 L ), including the "normal" LAEs. The 
luminosities are also in the same range as high redshift sub-mm 
galaxies (Chapman et al. 2005). At z = 0.3, roughly half of the 
galaxies lie in the range of normal star forming galaxies, seven 



Fig. 2. IR to Lya luminosities. Filled stars are sources at z — 
0.3, filled circles are at z = 2.3. The filled triangle is the only 
ULIRG detected in the shallow MIPS survey. Black indicates 
galaxies, and blue are unclassified LAEs. The solid diagonal line 
corresponds to a flux ratio of 0.02 %, other diagonal lines are 
0.05 and 0.005 %. Grey lines mark the two MIPS detection limits 
in the deep and wide surveys (vertical), as well as the flux limit 
in Lya for the z23 sample (horizontal). 

have LIRG luminosities (10 11 < Ljr < 1O 12 L ), and two have 
ULIRG luminosities. Note that at z — 2.3, detected sources are 
automatically ULIRGs, as the detection limit in the deep survey 
in the 24pm band corresponds to logL/R = 12.4 L , and in the 
shallow survey to logL/« = 12.9 L . 

In Fig.[2]the Lya luminosities are shown as a function of the 
infrared luminosities of the galaxies. Here and in the following 
analysis, we have chosen to be conservative and have excluded 
all AGN from the samples. In Sect. 14. 11 we return to the question 
of AGN and test how robust the results are against AGN inclu- 
sion. In Fig. |2] the z23 LAE candidates with MIPS detections 
seem to follow a given trend between the two flux measure- 
ments. The best-fit ratio between Lya and infrared luminosity 
is ~ 0.02 %. In the z03 sample, the Lya luminosity interestingly 
stays constant as a function of infrared luminosity. This indicates 
that the physical processes governing the Lya and the IR lumi- 
nosities are not related at low Lya and/or IR luminosities, al- 
though the relation seen at bright luminosities is based on small 
number statistics. As the blue points in this case (sources with 
no previous identification as either galaxy or AGN) are mixed in 
the population of galaxy LAEs based on Lya luminosity, these 
are hereafter considered as normal LAE candidates. 

The bolometric luminosity can be calculated from the in- 
frared and ultraviolet luminosity according to 

Lboi = Li R + Z4600, (2) 

where 

^1600 = ^160o/l60o4^rfL 2 . (3) 

This can further be converted to a dust unobscured star formation 
rate, assuming that the bolometric luminosity includes all the re- 
processed light from star forming regions (Kennicutt 1998): 

SFR = 4.5 x 10~ 44 x Lboi [erg s -1 ]. (4) 

The star formation rates found from the bolometric luminosity 
are in the range 500 - 5000 M yr" 1 for the non-AGN z23 LAEs 
and 2 - 3300 M yr 1 for the non-AGN z03 LAEs. Comparing 
the star formation rates found from the Lya line and from the 
bolometric luminosity, we find a median ratio of 0.0043+0.0025 
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Table 1. IR fluxes, star formation rates and derived dust attenuation. 



3 



z = 2.3 


LAE_ 




Lir 


Lboi 


S FRtoi 


A 1600 


Offset 


Remark 


COSMOS_# 




erg s" 1 


erg s 


M Q yr 1 




arcsec 




20 


154+16 


46.50 ± 0.05 


46.50 ± 0.07 


1420 ± 220 


5.51 ±0.18 


0.02 


AGN 


36 


835 ± 20 


47.69 ± 0.02 


47.69 + 0.02 


21830 ± 1000 


7.65 ± 0.05 


0.02 


AGN 


z = 0.3 


ID 


F^fm 


Lir 


Lbot 


S FRboi 


^1600 




Remark 






erg s~' 


erg s 


M Q yr" 1 








CDFS.1348 


1.63 ±0.03 


45.19 + 0.01 


45.20 ±0.01 


70.9 ±2.1 


4.91 ± 0.03 






CDFS.1821 


0.38 ± 0.05 


44.17 + 0.08 


44.20 ± 0.09 


7.1 ± 1.4 


2.69 ± 0.21 




Unclass. 



Columns are: (1) flux in the 8/24 yum data, (2) infrared luminosity, (3) bolometric luminosity, (4) star formation rate, (5) dust attenuation 
derived from the infrared luminosities, (6) offsets between MIPS catalogue source and LAE candidate in z23 sample, and (7) other remarks. 
The IDs for the z03 sample are from Deharveng et al. (2008). Full table is available in the online version. 



and 0.034 + 0.18 for the high and z03 non-AGN LAEs, where 
the error bars indicate the spread in the values. As the SFR found 
from the bolometric luminosity is the total SFR of the galaxy, 
tracing the same population of star forming objects as those cre- 
ating the Lya emission, the Lya escape fraction can be calcu- 
lated as the ratio between these two values. In the z23 sample, 
this escape fraction is ~ 0.4% for infrared selected Lya emit- 
ters, whereas it covers the whole range from a few to 100% in 
the z03 sample. Furthermore, for the z23 sample, five out of 
the 9 non-AGN LAEs have SFRs that should make them ob- 
servable in the Chandra X-ray images of the COSMOS field 
(logLx > 43.3 erg s ), if we assume the SFR-X-ray luminos- 
ity conversion of Ranalli et al. (2003). That they are not must 
indicate that the X-ray sources are obscured by dust. 

In Meurer et al. (1999), a relation between the ratio of in- 
frared to ultraviolet flux and the dust attenuation A1600 was de- 
rived: 

log = log(10°- 4xAl60 ° - 1) + 0.076. (5) 

^1600 

The spread in A 1600 for the galaxies here is large. The attenuation 
in both samples reaches as high as eight magnitudes in Ai6oo- In 
the z03 sample the attenuation spreads down to one magnitude, 
whereas the z23 galaxies all have Aieoo > 3. The distribution is 
flat in both samples. A plot of the SFR ratios against the dust 
attenuation derived according to this equation is found in Fig. [3] 
The SFR ratios, which are proxies of the Lya escape fraction, 
are seen to decrease with increasing dust attenuation, similar to 
what is seen in local galaxies (Atek et al. 2009). The appearance 
of this plot, however, stems from the relation between the x- 
and y-axis and Lya equivalent width (EW). Lines of constant 
EW have been overplotted in Fig. [3] All objects have EWs of 
roughly 20 - 400 A. 



4. ULIRG fraction and redshift evolution 

The finding that Lya emitting galaxies can also be very IR- 
bright, hence dusty, is not a novel result. Chapman et al. (2003, 
2005) showed that a set of sub-mm galaxies had Lya in emis- 
sion. Very red LAEs have also been found in other surveys. At 
Z ~ 3.1, Nilsson et al. (2007; N07, logL Lw > 41.88 erg s" 1 ) 
found one red LAE in a sample of 24 and Lai et al. (2008; 
L08, \ogLi ya > 42.07 erg s _1 ) presented four red LAEs among 
162, also at z ~ 3.1. To determine the fraction of ULIRGs in 
the z ~ 3.1 surveys, we studied the ULIRG template SEDs of 
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Fig. 3. Star formation rate ratios between Lya and infrared de- 
rived SFRs versus A^oo- Stars are sources at z = 0.3, and filled 
circles are sources at z — 2.3. The solid lines follow the rela- 
tion of constant Lya restframe equivalent width, with the lower, 
middle and upper lines representing EW - 20, 100, 400 A. 

Vega et al. (2008) and find that z = 3.1 ULIRGs must be de- 
tected in the Spitzer IRAC bands to the flux limits of ECDF- 
S. Furthermore, the colours in the IRAC bands have to be red; 
R - CM(8yum) > 2.5,C/zl(3.6yum) - Ch4 > 2. In the N07 
sample, only one object is detected in the IRAC bands, with 
R - Ch4 = 6.5, resulting in one ULIRG in this sample (1/24). 
In the L08 sample, 18 galaxies are detected in the IRAC bands. 
However, the colours of these objects, including the red objects, 
are all Chi - ChA < 2 (their Fig. 2). As a result, no ULIRGs 
are found in their survey (0/162). Combining these two samples, 
the number of ULIRGs at z ~ 3.1 are 1 out of 186, resulting 
in a ULIRG fraction of 0.5^ 5%, where the lcr error bars are 
based on Poisson statistics. Finally, Colbert et al. (2006) found 
that three out of 22 z ~ 2.3 LAEs (14^ 3 %, excluding their de- 
tections of ULIRGs in Lya blobs[3, see also Francis et al. 2001) 
were bright in the observed infrared, and that they were ULIRGs. 
The percentage of ULIRGs in the z = 2.3 sample presented here 
is 14^%, if only the objects in the deep survey are considered 
(7/50). This number is a lower estimate, as the sensitivity of the 
deep COSMOS MIPS data only reaches logL /K = 12.4 L at 

1 It is unclear what mechanisms power Lya blobs, and if they are a 
distinct population of objects or simply the tail of the size distribution 
of LAEs. Due to these uncertainties, we chose a conservative selection 
and disregard these sources here. 
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Fig. 4. ULIRG fraction in Lya samples as a function of age of the 
Universe. The points at z ~ 2.3 are lower limits since both sur- 
veys were incomplete in infrared luminosity. A transition from 
nearly no ULIRGs to ~ 10 - 20% ULIRGs in Lya selected sam- 
ples is seen at z ~ 2.5. 

z = 2.3. In the z = 0.3 sample, the results are complete in in- 
frared luminosity, and the percentage of ULIRGs in the sample 
is 20^g 2 % (6/30), where we again have included 6 unclassified 
objects but not the AGN. We return to the robustness also of this 
result against AGN inclusion in Sec. l4.1l below. 

In Fig. [4] the ULIRG fractions are plotted as a function of 
redshift. The data points are few, but there is a clear trend for the 
ULIRG fraction to grow from high redshifts up to the present 
day. There even seems to be a hint that the transition from al- 
most zero ULIRG fraction to the current value is rather sudden. 
For the sole purpose of illustrating the time and steepness of the 
transition, we overplot in Fig. [4] a hyperbolic tangents function 
of the form: 

UF(z) = ^(1 - tanh(0 (z - z, r ))). (6) 

In this equation, UF is the ULIRG fraction, 9 represents the 
steepness of the transition, and z lr is the transition redshift. The 
function plotted in Fig.|4]has UFq = 0.20 (present-day fraction), 
= 2.28, and z fr = 2.52. This means that the transition redshift 
for the ULIRG fraction is z ~ 2.5; the ULIRG fraction is negligi- 
ble at higher redshifts, whereas at lower redshifts it approaches 
~ 25%. 

4.1. Robustness of the result 

For the definition of the Lya-selected ULIRG subsamples we 
used the common definition based on L/«. IR-selected ULIRG 
samples are known to have a significant (15-50%) fraction of 
galaxies with AGN components, but it is also known that their 
total luminosity is often dominated by star formation rather than 
by the AGNs (Veilleux et al. 2009). Nevertheless, to make sure 
that our Lya selection did not bias the results, we chose the most 
conservative approach and therefore removed all X-ray detected 
objects from the samples. We now ask the question of whether a 
less conservative approach is possible, and if that would change 
the conclusions. We therefore searched for an additional test to 
certify that a candidate is a ULIRG. Using the ULIRG template 
SEDs of Vega et al. (2008), we find that a useful index at z = 2.3 
is optical R band minus Ch4 of IRAC. All Vega templates have 
R - Ch4 > 2.7, see Fig. [5] The points in the plot show all the 
LAE candidates at z — 2.3 from Nilsson et al. (2009) with > 3<x 
detections in R and Ch4. The black points are those detected with 
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Fig. 5. R - C/?4 colours for the LAEs at z = 2.3 with at least 
3cr detections in each band. For a description of the different 
data points, see text. The solid line marks our colour selection 
criteria (R - Ch4 > 2.7). 

MIPS, whereas the green points are not MlPS-detected. X-ray 
detected objects are marked with a pink ring. Four ULIRGs are 
not detected in Ch4 and are shown with upper limits. It can be 
seen from this plot that the flux selection with MIPS very cleanly 
samples the region of colour space where ULIRGs are expected 
to be, except the very brightest AGN that are too blue for this 
selection. Selecting ULIRGs based on the infrared flux and the 
R - Ch4 colour, and including all objects that have the correct 
colours within lcr, changes the z = 2.3 sample by increasing the 
number of ULIRGs from seven of a total 50 LAE candidates in 
the first selection to eight of a total 60 candidates in the second 
selection. The new fraction is 13*1%, in perfect agreement with 
the initial result. 

The ULIRG fraction at z — 0.3 was obtained using the 
Deharveng et al. (2008) sample but excluding the AGN reported 
by Cowie et al. (2009). If the AGN are included, the ULIRG frac- 
tion drops to 19^g 2 % rather than the previous 20*g 2 %, while ig- 
noring the Deharveng et al. (2008) sample. Using only the sam- 
ple of Cowie et al. (2009) gives a ULIRG fraction of 12 + | 6 %. All 
of those results are identical to within lcr errors, so the z = 0.3 
result is also robust against details of the sample selection. 

5. Conclusion 

We have here presented the infrared properties of two samples 
of Lya emitters, at z — 0.3 (24 objects) and z = 2.3 (25 objects). 
The samples were originally Lya-selected but here we consid- 
ered only the subsamples that were also detected with Spitzer. 
The two subsamples were selected to determine the fraction of 
powerful dusty starburst galaxies (ULIRGs) among LAEs and 
the redshift evolution of this fraction. Our most important con- 
clusions are the following. 

i) At all redshifts below z — 3, a non-zero fraction of LAEs 
are found to be ULIRGs. This result directly contradicts the clas- 
sic view that dust and Lya emission are mutually exclusive, and 
it holds together with the finding that Lya emission at high red- 
shifts correlate with metallicity (M0ller et al. 2004) makes the 
case that we must re-think the importance of dust for the Lya 
escape fraction. 

ii) There is evidence for a strong evolution in ULIRG frac- 
tion from redshift 3.1 to the present universe. The fractions de- 
rived have been shown to be very robust against different selec- 
tion criteria. This explains why there was severe disagreement 
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about the colour of Lya galaxies originally. At redshifts around 
three they were reported to be young, blue, low-dust starbursts 
(Warren & M0ller 1996, Fynbo et al. 2000), while at redshifts 
closer to two they were reported to be red/dusty (Stiavelli et 
al. 2001). 

Hi) There may be evidence that the evolution of the ULIRG 
fraction is not linear in time, but rather that there is a jump from 
redshift three to two, and then a more gradual evolution to the 
present time. Even though the evidence is merely suggestive at 
present, this jump does coincide with the number density evolu- 
tion of quasars and suggests a connection between those classes 
of objects. This question bears strongly on the connection be- 
tween starbursts and the formation of quasars, so it should be 
investigated more thoroughly via dedicated surveys of LAEs at 
several different redshifts in the range z = 2 to z - 3.5. We 
provide a simple formalism in the form of a three parameter 
function which is well suited to future tests of how sharp the 
transition is. 

From the accumulation of evidence, it is clear that some Lya 
emitting galaxies are very red and dusty, and vice versa that 
at least some ULIRG galaxies exhibit Lya emission. As was 
shown in e.g. Finkelstein et al. (2009), it is possible that for 
some LAEs the Lya emission is not affected by the dust at all, 
or at least not along the line-of-sight towards us. The possibility 
that Lya equivalent widths may be enhanced in dusty environ- 
ments has also been discussed in, e.g., Neufeld (1991), Hansen 
& Oh (2006), and Finkelstein et al. (2008, 2009). Interestingly, 
one of the ULIRG LAEs found at high redshift in the Nilsson 
et al. (2009) sample also has by far the largest equivalent width 
of their sample, with a restframe Lya equivalent width of nearly 
800 A. Considering that samples of LAEs now have for several 
years consistently included red, dusty galaxies, it is advisable to 
stop considering Lya emitters as one type of galaxy, as we have 
seen in this and previous publications (Nilsson et al. 2009) that 
Lya emitters exhibit a range of properties and that this range 
increases with decreasing redshift. 
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